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1. Introduction 

The nonet of light scalar mesons formed by a = /o(500), K = £q (800), «o(980) and / (980) 
is poorly understood. Compared to expectation all nine states are rather light and their ordering is 
inverted. For example in a standard quark antiquark picture the <3o(980) states, which have 7=1, 
must necessarily be composed of two light quarks, e.g. ao(980) = du, while the K states with 
I = 1/2 are made from a strange and a light quark, e.g. K = su. Consequently, K should be heavier 
than «o(980), since m s > m u ^. In experiments, however, the opposite is observed, i.e. m(fc) = 682± 
29MeV, while m(a (980)) = 980±20MeV %. On the other hand in a four-quark or tetraquark 
picture the quark content could be ao(980) = duss and K = su(uu + dd) naturally explaining the 
observed ordering. Moreover, certain decay channels, e.g. ao(980) — > K + K, indicate that besides 
the two light quarks also an ss pair is present and, therefore, also support a tetraquark interpretation. 
A detailed discussion of light scalar mesons can be found e.g. in [Q, H], 

In addition to the light scalar mesons there are also various tetraquark candidates among the 
heavy mesons, e.g. D* (2317) and 0^(2460) (cf. e.g. [|, |]) or the charmonium states X(3872), 
Z(4430), Z(4050), Z(4250). 

Here we report about the status of an ongoing long-term project with the aim to study possible 
tetraquark candidates from first principles using lattice QCD. Parts of this work have already been 
published [g^]. 



2. Lattice setup 

We use gauge link configurations with Nf = 2 + 1 + 1 dynamical quark flavors generated by 
the European Twisted Mass Collaboration (ETMC). 

We have studied several ensembles with the same rather fine lattice spacing a 0.086 fm. The 
ensembles differ in the volume (L/a) 3 x (T/a) = 20 3 x 48, . . . , 32 3 x 64 and the unphysically heavy 
light quark mass corresponding to m % 280 MeV . . .460MeV. Details regarding these gauge link 
configurations can be found in [g g 0, H 0. 



Currently we ignore disconnected diagrams, which are technically rather challenging (cf. the 
outlook in section ||). An important physical consequence is that the quark number and the anti- 
quark number are separately conserved for each flavor. Therefore, there is no mixing between uu, 
old and ss resulting in an r\ s meson with flavor structure ss instead of r] and X]' . 

Further lattice details and technicalities can be found in [^]. 



3. Four-quark creation operators 

In the following we focus on the ao(980) sector, which has quantum numbers I(J P ) = 1(0 + ). 
As usual in lattice QCD we extract the low lying spectrum in that sector by studying the asymptotic 
exponential behavior of Euclidean correlation functions 

C jk (t) = ((^-(OM(O)). (3.1) 

6 j and 0\ denote suitable creation operators, i.e. operators generating the #0(980) quantum num- 
bers, when applied to the vacuum state. 
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Assuming that the experimentally measured ao(980) with mass 980 ±20 MeV is a rather 
strongly bound four quark state, suitable creation operators to excite such a state are 



<fKK molecule 
«o(980) 

adiquark 

ao(980) 



= £^-(x) 75 «(x)j^(x) 75 ,(x)j (3.2) 
£ (e a V(x)Cy 5 <f ' r (x)) (e ade u d > T {x)CY 5 s e (x)^ . (3.3) 



The first operator has the spin/color structure of a.KK molecule (s(x)j5u(x) and d(x)j5s(x) corre- 
spond to a kaon K and an antikaon K at the same position x). The second resembles a bound diquark 
antidiquark pair, where spin coupling via C75 corresponds to the lightest diquarks/antidiquarks (cf. 
e.g. § 0,0). 

Further low lying states in this sector are the two particle states K + K and r] s + %. Suitable 
creation operators to resolve these states are 

^two-particle = ^, (x)75M(x) ^£J (y)75 , (y) ^ (3 .4) 
^particle = ^,( X ) 75 ,( X )^£J( y)75M ( y )). (3.5) 

4. Numerical results an their interpretation 

We first discuss numerical results for the ensemble with the smallest volume, (L/a) x (T/a) = 
20 3 x 48, which corresponds to a spatial extension of L 1.72fm. This ensemble is particularly 
suited to distinguish two-particle states with relative momentum from states with two particles at 
rest and from possibly existing £0(980) tetraquark states (two-particle states with relative momen- 
tum have a rather large energy because one quantum of momentum 
p min = 27r/L«720MeV). 

Figure [j]a shows effective mass plots from a 2 x 2 conelation matrix with a KK molecule 
operator ( ^2| ) and a diquark-antidiquark operator (|3~3|). The corresponding two plateaus are around 
HOOMeV and, therefore, consistent both with the expectation for possibly existing a , o(980) tetra- 
quark states and with two-particle K + K and r\ s + % states, where both particles are at rest 
(m(K + K) « 2m(K) 1198MeV; m(r\ s + n) m(ri s ) + m(n) w 1115MeV in our lattice setup). 

Increasing this correlation matrix to 4 x 4 by adding two-particle K + K and r\ s + % operators 
(eqs. (3^) and (|3~3|)) yields the effective mass results shown in Figure [j]b. Two additional states 



are observed, whose plateaus are around 1 500 MeV .. .2000 MeV. From this 4x4 analysis we 
conclude the following: 

• We do not observe a third low-lying state around 1100 MeV, even though we provide opera- 
tors, which are of tetraquark type as well as of two-particle type. This suggests that the two 
low-lying states are the expected two-particle K + K and r\ s + n states, while an additional 
stable ao(980) tetraquark state does not exist. 

• The effective masses of the two low-lying states are of much better quality in Figure |l|b 
than in Figure [j]a. We attribute this to the two-particle K + K and rj s + n operators, which 
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b) KK molecule , diquark , K+K two-particle , r\ s +n (4x4 matrix) 




c) operator content, state 



j = KK molecule 
j = diquark 
j = K+K two-particle 
j = n s +rc two-particle 



d) operator content, state 1 



j = KK molecule 
j = diquark 
j = K+K two-particle 
j = T] s +7i two-particle 



Figure 1: ao(980) sector, (L/a) 3 x (T/a) = 20 3 x 48. a) Effective masses as functions of the temporal 
separation, 2x2 correlation matrix (operators: KK molecule, diquark-antidiquark, eqs. (3.2) and (3.3)). 
b) 4 x 4 correlation matrix (operators: KK molecule, diquark-antidiquark, two-particle K + K, two-particle 
Tj v + n, eqs. (3.2) to (3.5)). c), d) Squared eigenvector components of the two low-lying states from b) as 
functions of the temporal separation. 



presumably create larger overlap to those states than the tetraquark operators. This in turn 
confirms the interpretation of the two observed low-lying states as two-particle states. 

• To investigate the overlap in a more quantitative way, we show the squared eigenvector com- 
ponents of the two low-lying states in Figure [T]: and Figure |l|d (cf. [ |13| ] for a more detailed 
discussion of such eigenvector components). Clearly, the lowest state is of rj s + n type, 
whereas the second lowest state is of K + K type. On the other hand, the two tetraquark 
operators are essentially irrelevant for resolving those states, i.e. they do not seem to con- 
tribute any important structure, which is not already present in the two-particle operators. 
These eigenvector plots give additional strong support of the above interpretation of the two 
observed low lying states as two-particle states. 

• The energy of two-particle excitations with one relative quantum of momentum can be esti- 
mated according to 

m(l+2,p = p min ) « ^m(l) 2 +/4 n + ^m(2) 2 +/4 n , p min = ^. (4.1) 

Inserting the meson masses corresponding to our lattice setup, m{K) « 599 MeV, 
m{r\ s ) « 774 MeV and m{%) 341 MeV, yields m(K+K,p = p min ) w 1873 MeV and 
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m{r\ s + 71, p = p m i n ) ~ 1853 MeV. These numbers are consistent with the effective mass 
plateaus of the second and third excitation in Figure [j]b. Consequently, we also interpret 
them as two-particle states. 

We obtained qualitatively identical results, when varying the light quark mass and the space- 
time volume, as discussed in section ||. Corresponding plots are shown in [[7|]. 

Using exactly the same techniques, i.e. four-quark operators of tetraquark and of two-particle 
type, we also studied the K"-sector (for details cf. [^]). Again we find no sign of any four-quark 
bound state besides the expected two-particle spectrum (in this case K + % states). Note that this 
result is in contradiction to a very similar recent lattice study of the K meson [16], where an addi- 
tional low lying four-quark bound state has been observed. 



5. Inclusion of singly disconnected diagrams 

As mentioned in section || disconnected diagrams have been ignored for the results presented 
so far. In this section we briefly discuss our strategy for computing such diagrams for specific 
four-quark operators. 

For correlation functions of four-quark operators with flavor structure qxqxq^qi, a so-called 
singly-disconnected diagram has to be computed. Tetraquark candidates with this flavor structure 
are e.g. the previously discussed £20(980) = ssdu (cf. Figure ||) or 
D; (2317),D s i(2460) = {uu + dd)cs. 

While for connected four-quark diagrams (i.e. all four quark propagators connect the timeslices 
at time and time t) standard point-to-all propagators can be used, applying exclusively such 
propagators is not possible in practice for singly disconnected diagrams. The reason is that one has 
to include a sum over space, £ x , at least on one of the two timeslices (wlog. at time t in Figure |2|), 
to project to zero momentum. This in turn requires a all-to-all propagator of quark flavor q\ on 
that timeslice, due to £ x <?i {x,t)q\ (x,t) . . . (the s quark propagator at timeslice t represented by the 
solid red line in Figure |||). 

Since all-to-all propagators are prohibitively expensive to compute, they are typically esti- 
mated stochastically. While using a single stochastic propagator for a specific diagram typically 
results in a favorable or at least acceptable signal-to-noise ratio, using a larger number of such 
propagators drastically increases statistical errors. Therefore, we decided for the following strat- 
egy: three quark propagators (the q\ -loop at timeslice and the q2 and 173 propagators connecting 
the timeslices and t) are realized by exact point-to-all propagators, while the remaining propaga- 
tor (the <7i-loop at timeslice t) is estimated stochastically, using random Z2 x Z2 timeslice sources. 

First numerical tests of this strategy performed for both £?o(980) and D* (2317) have been 
promising in the sense that the statistical errors of the singly disconnected diagrams for the molecule 
type operator ( |3.2| ) are of the same order of magnitude as the statistical errors of the corresponding 
connected diagrams, when investing a comparable amount of HPC resources. 



6. Conclusions and future plans 

We have studied the £?o(980) and the K channel by means of N/ = 2 + 1 + 1 flavor lattice 
QCD using four-quark operators of molecule, diquark and two-particle type. Besides the expected 
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two-particle spectrum (two essentially non-interacting pseudoscalar mesons) no indication of any 
additional low lying state, in particular no sign of a four-quark bound state could be observed. This 
suggests that both the a$ (980) and K meson have either no sizeable tetraquark component or they 
are rather weakly bound unstable states. 

To investigate the latter one needs to study the volume dependence of the two-particle spectrum 
in the corresponding sectors ("Luscher's method", cf. e.g. [ ]T7[ |D|,[l9|]). Such computations are very 
challenging using lattice QCD, but first results have recently been published (cf. [2^, 21]). We plan 
to perform similar computations with our setup in the near future. 

Moreover, certain possibly present systematic errors need to be studied, quantified and re- 
moved: (1) disconnected diagrams have to be computed and included (cf. section |J); (2) lattice 
discretization errors and the continuum limit has to be studied; (3) computations at even lighter 
and, therefore, more realistic u/d quark masses would be desirable. 
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